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1 Introduction. 



The adequate description of the gravitation field of an astrophysical object 
has been an important subject in both relativistic and Newtonian gravity 
since their origins. The particular case of the gravity associated to axially 
symmetric bodies has played a central role in this discussion. Recently, Mer- 
rit rjTJ| found, from detailed modeling of triaxial galaxies, that most of the 
galaxies must be nearly axisymmetric, either prolate or oblate. In Newto- 
nian theory the gravitational potential of axially symmetric bodies can be 
always represented by its usual expansion in terms of Legendre polynomials 
(zonal harmonics) . In general relativity we have that the solution of the vac- 
uum Einstein equations associated to static axially symmetric bodies has a 
simple form with only two metric functions 0. 

As a consequence of the self interaction of the general relativistic grav- 
itational field we have a rather elaborate description of multipolar fields, the 
covariant multipole moments ||. For the special case of axially symmet- 
ric spacetimes we have geometric properties that can also be used to define 
a class of multipole moments [|, It is highly non trivial to prove the 
equivalence between these two definitions ||. Also, due to the nonlinearity 
of the Einstein equations we do not have a simple expansion of the met- 
ric in terms of relativistic multipoles like the expansion of the potential in 
Newtonian gravity. Moreover, it is not know how to construct a solution to 
the vacuum Einstein equation with assigned relativistic multipoles. Also the 
physical consequences of having different relativistic moments are rather ob- 
scure. Furthermore, the problem to model a static axially symmetric center 
of attraction in general relativity with some multipolar perturbations that 
have a prescribed Newtonian limit is not a well define problem, as a matter of 
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fact we have an infinite number of solutions of the vacuum Einstein equations 
with the same Newtonian limit. This problem arises in a natural form when 
one studies the orbits in deformed centers of attraction, see for instance 0. 

In this paper, in order to shed some light in the physical meaning of 
the relativistic multipolar expansions, we consider different static solutions of 
the axially symmetric vacuum Einstein equations that in the non relativistic 
limit have same Newtonian moments. The particular solutions considered 
represent a large class of rather natural solutions that are generalizations of 
some well know ones. 

Along this article we shall use the terminology "distorted black hole" 
introduced in || to indicate an attraction center with multipolar moments, 
examples are: a) A true black hole (or a dense object) surrounded by a distri- 
bution of matter like a ring or a small disk formed by counterrotating matter, 
i.e., built by approximately the same number of particles moving clockwise 
as counterclockwise. Even though, this interpretation can be seen as a device 
to have a static stable configuration there are observational evidence of disks 
made of streams of rotating and counterrotating matter ||. b) An axially 
symmetric static dense object with either polar deformations or polar jets. 
In our analysis the horizon will play no role, therefore our results will be 
valid for any nearly static axially symmetric attraction center. 

Since the advance of the perihelion is one of the most significant fea- 
tures of the general relativistic celestial mechanics will be considered in some 
detail. We find discrepancies in the fourth order of the dimensionless pa- 
rameter (mass of the attraction center) /(semilatus rectum). We present an 
evolution equation for the difference of the radial coordinate due to the use 
of different general relativistic multipole expansions. 

In Sec. 2 we present a summary of the main expressions associated to 
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the Weyl solutions. In Sec. 3 we study the Newtonian limit of four different 
solutions to the Einstein equations. Also we study the relativistic multipoles 
moments for each solution. These solutions are the Erez-Rosen-Quevedo 



(ERQ) solution l lOfl , a solution characterized by the usual Newtonian mul- 
tipoles and generalizations of both. We present two graphics that show the 
difference between the solutions. In the next section, Sec. 4, we study series 
expansions of the solutions restricted to the equatorial plane. Also several 
differences of the metric potentials are considered. In Sec. 5 an evolution 
equation for the difference of the radial coordinate due to the use of different 
multipole expansions is presented. We solve this equation for circular orbits 
in the equatorial plane. In Sec. 6, the perihelion shift is computed to the 
order in which the first discrepancies appears due to the use of the different 
multipolar expansions. As a particular case we have the perihelion shift in 
a Schwarzschild spacetime up to a fourth order. We also solve the difference 
equation of the preceding section for generic orbits on the equatorial plane. 
Finally, in Sec. 7, we present a short discussion of our main results. 

2 Static axially symmetric solutions of Ein- 
stein vacuum equations. 

The external gravitational field of produced by an axially symmetric body 
can be well described by the Weyl metric 0, 

ds 2 = e 2 ^dt 2 - e- 2 ^{r 2 dv 2 + e 2 \dr 2 + dz% (1) 

where the functions ip an d 7 depend only on r and z; the ranges of the 
coordinates (r, ip, z) are the usual for cylindrical coordinates and — oo < t < 
+oo. The vacuum Einstein equations {R^ u = 0) reduces to the usual Laplace 
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equation in cylindrical coordinates, 

^,rr + 1p,r/r + ^,zz = 0, (2) 

and the quadrature, 

dyfy) = r[(f r - f z )dr + 2^^ iZ dz\. (3) 

When ip satisfies the Laplace equation this differential is exact. 

The Schwarzschild solution in Weyl coordinates takes the form 

1, R+ + R. -2m 
ib s = -In — , (4) 

^ b 2 R + + +2m' V ; 

1 (R + + R^) 2 -4m 2 
75 = 2 M R-Ri ' (5) 

where 



R ± = ^r 2 + (z ± m) 2 . (6) 

The function ips is just the Newtonian potential of a bar of length 2m and 
density 1/2. The fact that in Weyl coordinates the Schwarzschild metric 
does not look spherically symmetric and that the horizon is squeezed into a 
line of length 2m has prevented the wide use of these coordinates. Since the 
superposition of static axially symmetric perturbation can be performed in a 
simple way and exact results can be obtained, we believe that is worth to use 
these coordinates. More suitable coordinates are the spheroidal coordinates, 

x = (R+ + R_)/(2m), y = {R + -RJ)/{2m). (7) 

Note that the coordinates (x, y) are dimensionless and 1 < x < oo and 
— 1 < 2/ < 1. The coordinate x is essentially a radial coordinate and y the 
cosine of an angle. The inverse of (0) are, 



m 



x 2 — l)(l — y 2 ), z = mxy. (8) 



In this case the Weyl metric takes the form 

dx 2 dy 2 
"x 2 — 1 1 — y 



ds 2 = e 2 ^dt 2 - m 2 e- 2 ^(x 2 - y 2 ){^- + -^~ 2 ) + (x 2 - 1)(1 - y 2 )d^% 



(9) 

and the Einstein equations, 



T,a: = ^TZ^i^ 2 - - x(l - y 2 )i/>l - 2y(x 2 - l)iPM (ii) 

t,, = irri^^ 2 - x )€ - rt 1 - f>i + 2y(i - iflVw^l- (12) 



[(i-^M a »],x + [(i/ 2 -iM tf ],» = o J (10) 

a; z — 
x 2 -l 
x 2 — y 2 

In this coordinates the Schwarzschild solutions takes the particularly simple 
form 

^5 = ~Qo(x), (13) 

1, x- 1 

= -In , 

2 x + 1' 

1 x 2 — 1 

7. = 2 ln ^7' (14) 

From the coordinate transformations, a; = R s /m — 1, y = cos'i?, and (^) 
with (^) and flliD we recover the Schwarzschild metric in its usual form, 

ds 2 = (1 - 2m/R s )dt 2 - (1 - 2m/R s )- 1 dR 2 s - i? 2 (A? 2 + sintf 2 ^ 2 )- (15) 

We shall consider solutions of the form, 

ip = ip s + ip, (16) 

where ?p(x,y) represents the superposition of external multipolar fields solu- 
tions of Laplace equation with no monopolar term. The function 7 in this 
case can be written as 

7^]= 7s + 7, (I 7 ) 
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with 



7 = 27^5,^+7^], (18) 
dl[tps, i>] = r[{ips,rl,r ~ ips,zl,z)dr + {lp S ,rl,z + ips,zl,r)dz). (19) 

In this case the Weyl metric takes the form, 
rfs 2 = (i_^)e^ dt 2_f!!^ (20) 

that has been interpreted as a "black hole" with multipolar deformations [§]. 

Note that this metric for a given function ip - solution of the usual 
axially symmetric Laplace equation - and its associated potential 7 given by 
Eq. ([TJ) is an exact solution of the Einstein vacuum field equations. We 
shall consider several different classes of approximations. In principle, since 
we start with an exact framework we will be able to control the validity of 
the approximations. The relation of this Weyl approach in the case linear 



perturbations with the Regge- Wheeler [11 1 formalism - general linear per- 
turbations of a black hole - can be found in In the present paper, in 
some cases, we will go beyond the linear case. 



3 Multipoles and multipolar fields. 

In order to gain some understanding of the solutions of the vacuum Einstein 
equations with different relativistic moments, in this section, we study dif- 
ferent Weyl solutions that represent static attraction centers that have the 
same Newtonian limit. 

Following Ehlers ||14|| , for a given Weyl solution we define its Newto- 
nian limit as liniA=o ^(A, r, z)/\, where A = c~ 2 and c the light velocity, e.g. 



7 



for the Schwarzschild solution we put m = XGM in (^) and with the aid of 
l'Hopital rule we get 

limVs(A,r,*)/A = —GM/R. (21) 

A — 

with 

R = Vr 2 + z 2 . (22) 

There are an infinite number of different Weyl solutions that represent 
fields of multipoles with the same Newtonian multipolar limit. We shall 
restrict our study to two known solutions that have been used recently in 
some applications and other two that are generalizations of the formers. They 
are: 

a) The usual Newtonian multipolar fields, 

V> =E<fc Rk+ i » ( 23 ) 

k=l 

we have the same ip N as its Newtonian limit. Note that qt scales with A as 

b) Separating variables in (0) one finds the ERQ solution |15|, |10|1 , 



i- ERQ -t^0r^Q^) p ^v), (24) 

where Qk{x) and Pk{y) are the Legendre functions of the second kind and 
the Legendre polynomials, respectively. We have chosen the coefficients in 
the series to have 

hm^ Q (A,r,z)/A = ^, (25) 

where 4>n R ® an d are the n-term in the series (|24|) and (|23|), respectively. 

To find the limit is useful the Letelier identity ||16|| , 

^ / x / % 1 f m P n (a/m) , 

Qn(x)P n (y) = - / ^ I > da. (26) 

l J-m J r 2 + [z — a) 2 



c) The solutions of Laplace equation: 



u f r P n ((z + am)/R +a ) ( P n ((z - am) / R_ a ) 
i> =2^<lk[a + (l-a) j, (27) 

k=l rL +a ^-a 

00 (Oh + I V 

^%?24^^ W?fcW ' (28) 



x a = (R +a + R_ a )/(2ma), y a = (R +a - R_ a )/(2ma), (29) 



where 



R +a = yj r 2 + ( z + am) 2 , i?_ a = ^r 2 + (z - am) 2 , (30) 

a and a are arbitrary constants that we shall take as positive and less than 
one, we shall comeback to this point later, ip 1 is obviously a solution of the 
Laplace equation. To prove that ip 11 is also a solution of this last equations 
one can use the identity (26). -ip 1 and -ip 11 are nontrivial deformations of the 



Newtonian multipoles and the ERQ solution, respectively. Note that in the 
superposition (|16| ) the Schwarzschild term ips is not changed. Letting a = 
in ip 1 we get ip N > a l so putting a = 1 in ip 11 we recover ip ER ®. The global 
properties of these new solutions will be presented elsewhere. 

The first two solutions are the axially symmetric multipolar expansion 
for the Newtonian Laplace equation in spherical coordinates (usual expan- 
sion) and in spheroidal prolate coordinates, respectively. The second one 
has the property of having the monopolar term proportional to the Weyl 
potential that gives rise to the Schwarzschild solution [cf. Eq. (|T3"D]. The 
Schwarzschild solution with multipolar deformations represented by either 
expansions have been study by several authors (see for instance, JO], |17|1), 



specially in the important case of quadrupolar deformations [§, 15 . For re- 



cent applications see M, ff8l. The last two multipolar expansions, ip 1 and ip 11 , 



are new ones and are closed related to the previous ones. They have some 
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new parameters that do not appear in the Newtonian limit. They will be 
considered in order to explore the possibility to have relativistic multipoles 
not completely determined by its Newtonian ones. We shall comeback to 
this point later. 

The relativistic multipoles m k for a Weyl solution can be computed 
from the function ip evaluated on the axis of symmetry, in the following way 
[ PH , [T3[l : First we define the function 



£{x,y)= 1 + e2#E|y) (31) 

and find its value on the axis of symmetry r = 0, i.e., y — 1; so x — z/m. 
Then we calculate the quantities 

mk = (k + l)\ dz^ ls =° (32) 
with z = 1/z. The relativistic moments are given by 

m k =rh k + d k , (33) 

where the first six constants d k are fl3| : 



d = d\ = d 2 = d 3 = 0, 
c?4 = fh (rhl — fh 2 fh )/4:, 

d 5 = rh (m2Thi — m 3 mo)/3 + {fhi{fn\ — m 2 ?^o)/21. (34) 

The direct computations of the relativistic multipolar moments up to 
the sixteenth pole gives: 

a) For the Newtonian multipolar field ip N , 

m = m, m 1 = q h m 2 = q 2 , m 3 = g 3 - m 2 q 1 , (35) 
m 4 = g 4 — (8m 2 q 2 + 6mql)/7. 
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b) For the ERQ multipolar field i/> ERQ , 

m = m, m 1 = q 1 , m 2 = q 2 , m 3 = q 3 - 2m 2 q 1 /5, (36) 
m 4 = q 4 — 2{m 2 q 2 + 3mg 2 ) /7. 

c) For the multipolar field t[) T , 

m = m, m 4 = q±, m 2 = q 2 + 2(1 — 2a)amq±, (37) 
m 3 = ?3 + (3a 2 — l)m 2 gi + 3a(l — 2a)mq 2 
m 4 = q 4 + 4a(7a 2 - 4)(1 - 2a)m 3 q 1 /7 

+4(1 - 2a)amq 3 + 2(21a 2 - 4)m 2 g 2 /7 - 6mg 2 /7. 

d) For the multipolar field tjj 11 

m = m, m 1 =q 1 , m 2 = q 2 , m 3 = q 3 - (1 - 3a 2 /5)m 2 q 1 , (38) 
m 4 = q 4 — 6mqf/7 — 2(4 — 3a 2 )m 2 q 2 /7. 

Since our main goal is to study the physical effect of having different 
relativistic multipole moments is instructive to compute the differences of 
these multipoles for the above mentioned, solutions. We shall compute the 
quantities, 

Am k = m k - m%. (39) 

We find for the ERQ, I and II solutions Am = Ami = 0, and 

a) For the ERQ solution: 

Am 2 = 0, Am 3 = 3mgi/5, Am 4 = 6m 2 q 2 /7. (40) 

b) For the solution I: 

Am 2 = 2a(l — 2a)mqi, Am 3 = 3a 2 m 2 q\ + 3a(l — 2a)mq 2 , 

Am 4 = 4a(l - 2a)(7a 2 - 4)m 3 q 1 /7 + 6a 2 m 2 q 2 + 4a(l - 2a)mq 3 .(Al) 
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c) For the solution II: 



Am 2 = 0, Am 3 = 3a 2 m 2 qi/5, Am 4 = 6a 2 m 2 g 2 /7. 



(42) 



Along all the paper we shall compare the different solutions with ip N , of 
course this choice is arbitrary. 

Since the quadrupolar strength q 2 plays an important role in the above 
differences - moreover when qi = - we shall examine closely the quadrupo- 
lar potentials. A significant quantity is the quadrupolar deviation, 



that we shall study graphically. 

In Fig. 1 we present the deviation of tp 2 = i>l on the plane d = it/ 2 
for a = 1, and different values of the constant a: 0.25 (top curve), 0.5, 0.75, 
1.0 (bottom). We can have large deviations that asymptotically vanish. In 
Fig. 2 we show the deviation ip2 = V^ 7 011 ^he plane $ = ir/2 for different 
values of the parameter a: 0.25 (top curve), 0.5, 0.75, 1.0 (bottom). The 
case a = 1 corresponds to the ERQ solution. It is not difficult to show using 
Letelier identity that in the limit a = 0, ip 11 = ifj N , in other words the ERQ 
solution can continuously deformed in to ip N , fact that is confirmed in Fig. 
2 for the quadrupolar term. 

4 Solutions and differences. 

For our purposes, it is enough to consider the the multipolar series only up to 
the fourth order term, ip^. Since one of our objectives is the study of orbits 
of test particles moving around a deformed center of attraction, the inverse 
"radial coordinate", u = Rg 1 = [m(l +a;)] -1 , will be useful. Moreover, all 



(^2-OM 



(43) 
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the metric functions will be computed in the plane d = tc /2 up to the sixth 
order in u; this order will be the adequated for consistence in the expansions. 
Also we shall disregard the dipolar moment, q\ = (we shall comeback to 
this point latter). We find for the different functions ^ = ip 2 + ^3 + ^a- 

13 3 
ip N = --q 2 u 3 - -mq 2 u 4 - -(I0m 2 q 2 - q^u 5 
Z Z o 

5 

— m(lAm 2 q 2 - 3g 4 )u 6 , (44) 
8 

13 3 

^ERQ = _ u 3 _ m(?2M 4 _ (Q Am 2 q2 _ 7q 4 ) u 5 
Z Z 00 

_ JL m ( 80m 2 g 2 - 21q 4 )u e , (45) 
56 

13 3 
W = ~-q 2 u 3 - -mq 2 u 4 - -[2(5 - 3a 2 )m 2 q 2 + 4(2a - l)aq 3 - q 4 ]u 5 
Z Z 8 

5 

--m[2(7 - 9a 2 )m 2 q 2 + 60(2a - l)amq 3 - 3g 4 ]w 6 , (46) 
8 

13 3 
ip H = --Q2U 3 - -mq 2 u A - —[(70 - 6a 2 )m 2 q 2 - 7g 4 ]w 5 
Z Z 5o 

-— m[2(49 - 9a 2 )m 2 q 2 - 21g 4 )u 6 . (47) 
56 

The respective functions 7 are: 

7^ = -^u 4 - 3m 2 q 2 u 5 - -{3q 2 + 70m 3 q 2 - 5mg 4 )w 6 , (48) 
4 8 

*erq = _3^4 _ 3m 2 gaM 5 _ (21^ + 460m 3 g 2 - 35mg 4 )w 6 , (49) 
4 56 

3 1 

7 7 = --u 4 - 3m 2 q 2 u 5 - -{3q 2 2 + 70m 3 q 2 - 30a 2 m 3 q 2 - 5mg 4 )u 6 ,(50) 

4 8 

3 1 

7 7/ = — m 4 - 3m 2 q 2 u 5 (21q 2 + A90m 3 q 2 - 30a 2 m 3 q 2 

4 56 

-35mg 4 )w 6 . (51) 

The exact function 7 for the Schwarzschild solution with Newtonian multi- 



poles can be found explicitly for the general case |I7H, and in the case of the 



multipoles ijj ER Q in ||10| . These relations are rather formidable and of not 
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much use. An integral representation of the function 7 for the ERQ solution 
can be found in [IB]. For the solution / we were able to find 7 in the general 



case. For the solution // up to the quadrupolar moment. These results are 
rather cumbersome an will be presented elsewhere. 

Note that to eliminate conic singularities along the axis of symmetry 
we need to impose the boundary condition lim r= ^ = outside the source. 
The solutions with this boundary condition are asymptotically flat. 

Since our main goal is to study different solutions with the same New- 
tonian multipole moments we shall compute the quantities, 

5ip = -0 - ip N , 57 = 7-7^. (52) 
Note that Sip = 5ip and £7 = 5j. We find, 

Stf} ERQ = —m 2 q 2 u 5 + —m 3 q 2 u e , (53) 
28 28 

3 

dip 1 = -a[3amq 2 — 2(2a — l)q 3 ]mu 5 
15 

+— a[3amq 2 - 2(2a - l)g 3 ]mV, (54) 
Sip 11 = a 2 5iP ERQ . (55) 

And for the function 7, 

5>f«Q = V, (56) 
6^ = ^a[3amq 2 - 2(2a - l)g 3 ]mV, (57) 

fry**? = ^a 2 m 3 q 2 u 6 . (58) 
28 

These quantities will play and essential role in our analysis, we note that the 
quadrupolar strength is the main parameter that appears in these differences. 
Also in dtp 1 and 5j we have the octopolar constant, ^3. Note that the 
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solution if) 1 is highly asymmetric due to the weight constant a. For particles 
orbiting a central body deformed by this highly asymmetric multipolar field 
we might have instabilities, for simplicity we shall take the symmetric case, 
a = 1/2, in this case, 

9 45 
Sip 1 = -a 2 mq 2 u 5 + — 3a 3 m 3 q 2 u 6 , (59) 

Sy = ^- a 2 m 3 q 2 u 6 . (60) 

In summary, we have that the discrepancies in the solutions appear in the 
fifth order of the inverse radial variable u for the potential ip and the sixth 
order for the potential 7, 

i = ~ u 3 - \mq 2 if + (-(3m 2 q 2 + ^g 4 )u 5 + 0(u 6 ), (61) 
2 2 8 

7 = -^mq 2 u 4 -3m 2 g 2 w 5 + 0(u 6 ). (62) 

with 

P N = ^ , P ERQ = y , P 1 = |(5 - 3a 2 ) , = |(70 - 6a 2 ). (63) 
For the differences up to the fifth order in u we have, 07 = £7 = and 

Sip = 5ip = -5(3m 2 q 2 u 5 } (64) 

with 

5 pERQ = _jL 5P J = --a 2 , 5 p n = --a 2 . (65) 
H 28' P 4 ' H 28 V ; 

For general linear perturbations (not necessarily static axially sym- 
metric) of the Schwarschild black hole we have the Regge- Wheeler [|ll|] formal- 



ism. One can show ]12| that the simple case of linear static axially symmetric 



perturbations in Weyl coordinates can also be put in the Regge- Wheeler form, 
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as long as, the solutions do not have conic singularities. In the framework of 
above mentioned formalism one can also show that the dipolar perturbations 
can be eliminated by a transformation of coordinates as in the usual New- 
tonian case. This is not a trivial fact since relativistic dipolar moments are 
invariant tensorial quantities and as such cannot be made zero by a coordi- 
nate transformation. In the linear approximation the dipolar perturbations 
turn to be gauge dependent. 

5 An evolution equation for differences. 

A first and simple approach to the problem of looking for the discrepancies 
in the orbits of test particles moving around a deformed attraction center 
modeled with different relativistic gravitational fields (metrics) with the same 
Newtonian limit, is to compare circular orbits by means of an equation for 
the differences of the radial coordinate. The geodesic equation for the Weyl 
metric has three constants of motion, that in the plane y = determine 
completely the orbit of test particles: 

x 2 y 2 



e ^2_ m 2 e 2 7 -2^2_ y 2 )( _^_ 

x z — 1 i — y 
-mV 2 % 2 -1)(1- y 2 )v 2 = 1, (66) 

e 2 ^i = /, (67) 

e- 2 ^x 2 -l)(l- y 2 ) l f = — 2 , (68) 

where, as usual, i = dt/ds, etc., / and h are constants related to the test 
particle energy and angular momentum, respectively. 

From (|66|) -(68) we find the equation of the orbits in the equatorial 
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plane, y = 0, 

L2 2 ^2 2^ 

/ 2 e- 2 ^ - e 2 ^- 2 -^—x* - IL-I— = 1, (69) 
m z [X 2 — l) d m z x l — 1 

with x' = J. It is more convenient to write this equation in terms of the 
potentials, ip an d 7, and the variable u = Rg 1 = [m(l + x)] -1 . Hence, 

u « + u 2 e -2^_^ e -2(^) + J, e -2(7+^)_^ e -2(^)_ 2mM 3 e -27 = Q? (7Q) 

that has the general form, 

m /2 = F(m, 7 (m),^(u)), (71) 

where 

F = e~ 2 H-u 2 + 2mv? + ^e~^ - ie" 8 * + ^wT 8 *). (72) 

AT /l^ /z/ 



By differentiation of ( ffl] ) we obtain the general equation for the difference of 
the inverse radial coordinate of equatorial orbits, 

<5w'-— — 5m=— — 5 7 + — r<^h 73 
v 7 2m' 9m 2m' \ d*f dip ) 

where £7 and Sip represent the differences defined in (|52"|). 
Note that 

% = -I?- («) 
For circular orbits, u = uo, (i?s = 1/mo, the constant radius), we have 
F = 0. From fl7|) and (|73|) we get 

S^-(^) * ,75) 



dip du 



<<() 
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First, we shall compute the leading term in u of this last equation, from (|64|) 
we find, 

5u = ±(l-2l 2 )q 2 45(3. (76) 

This last equation in terms of the Schwarzschild like radial coordinate 
Rs reads, 

-SRs = l(l-2l 2 )q 2 ^. (77) 

Then the difference decays with the square of the radius. It is also instructive 
to express equation ( ff6| ) in terms of the dimensionless parameter, 

half Schwarzschild radius 



orbit radius 

We get, 



(78) 



5u = 5/3(1 -2l 2 )^u e 3 . (79) 
am 6 

For planets like Earth, Mercury, and Mars e is a little parameter e ~ 10 -8 . 
But, for a small test body orbiting around a neutron star we can have e ~ 
10 _1 . In this last case the multipolar fields can be originated, for instance, 
by polar jets of matter ejected by the neutron star. Finally, we want to 
comment that Eq. ( [75|) gives a simple expression involving the exact metric. 

6 The perihelion shift. 

To compare the different multipole expansions we shall study the perihelion 

shift of a test particle orbiting in the equatorial plane of the deformed "black 

hole". In the order of approximation used in the present work we can put 

= 1 + if) and e 1 = 1 + 7 in flTOj ) with no error, 

/2 2 I 2 — 1 2mu o 

u + u — — 2mu = 

n z n z 

( 9 I 2 - 1 2mu A 7 / 2l 2 1 2mu\ , , 

2 7 [u 2 — - 2mu 3 +2tf)i -— + — — — — . (80) 
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Doing ip = 7 = in the previous equation we obtain the well known equation 
for the motion of a test particle in the Schwarzschild spacetime. 

Since the first term in the series expansion of the functions 7 and ip 
are proportional to u A and u 3 , respectively [cf. Eqs. ( ^) - (|51]) 1 ) to use the 
orbital equation fl30|) without error we must limit the expansion of 7 and Tp 
to the order seven and five in the variable u, respectively. Therefore we shall 
consider all the series up to the fifth order in u, that is the order wherein the 
first discrepancies between the different multipole expansions appear. From 
Eq. (|80| ) we obtain, 



u' 2 + u 2 = [I 2 - 1 + 2mu + 2mh 2 u 3 + -qm 3 u 3 + -rm^u 4 + -sm 5 u 5 }/h 2 , (81) 

o A o 

where, 

« = \ ( 2< ' 2 5- <82) 



r 



(l5(/ 2 -l) + 10)^, (83) 
((4/3 + Q)(l 2 ~ 1) + 2/? + 9) *| - -(2l 2 - 1) * 



5 



4) 



2 

Note that q, r, and s are dimensionless quantities. The semilatus rectum for 
the usual Keplerian orbit of a particle attracted by a central body with no 
multipolar deformations is, p = h 2 /m. We shall use a dimensionless param- 
eter similar to the homonymous one of the preceding section, 

m 2 m half Schwarzschild radius 
h 2 p semilatus rectum 

Note that for a circular orbit with u = uq we have uq = 1/p. It is enlightening 
to write the equation (|3lD in terms of the new dimensionless variable w = pu, 

I 2 — 1 2 1 2 

w ' 2 + w 2 = p 2 —— + 2w + 2ew 3 + -qe 2 w 3 + -re 3 w A + -se 4 w 5 . (86) 
h l 3 2 5 
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By derivation we obtain, 

w" + w = 1 + 3ew 2 + qe 2 w 2 + re 3 w 3 + se 4 «A (87) 

We remind that the discrepancies between the different multipole expansions 
are controlled by the parameter j3 that it only appears in s (last term). We 
shall search a solution of equation ( |8~TD in the form of a series expansion up 
to the fourth order in e, 

w = wo + w\e + w 2 e 2 + w 3 e 3 + u>4e 4 , (88) 

From (|S7|) we find the system of linear equations, 

Wq + w = 1, 

w" + wi = 3wq, 

w' 2 ' + w 2 = Qw wi + qwl, 

w'l + w 3 = 3(2w w 2 + w\) + 2qw wi + rw^, 

w'l + w 4 = 6(w w 3 + W1W2) + q(2w w 2 + w\) + 3rwlwi + sWq. (89) 

The solution of this system can be expressed in terms of elementary functions, 
we find, 

wo = 1 + e cos <p, 

wi = 3(1 + e 2 /2) + 3e(psmtp — (e 2 /2) cos2(^, 

w 2 = 9(2 + e 2 ) + q{\ + e 2 /2) - (9/2)e<^ 2 cos (p - (3e 2 + ge 2 /6) cos 2<p 
+ (3e 3 /16) cos3</? + (27e/2 + 15e 3 /4 + qe)<p sin if — 3e 2 uj sin 2ip, 

w 3 = 135 + 81e 2 + 57e 4 /8 + 12q + 6e 2 g + r + 3e 2 r/2 - (81e^ 2 /2 

+45eV /4 + 3eg<^ 2 ) cos p - (27e 2 + 59e 4 /16 + 2e 2 g + e 2 r/2 

-9e 2 <p 2 ) cos2</? + (9e 3 /4 + e 3 g/8 - e 3 r/32) cos3</? - ^-e 4 cos4</? + 

16 
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(189eyV2 + 135eV/4 + 9eg</? + 5e 3 q<p/2 + 3ery>/2 + 3e 3 rp/8 

-9e<p 3 /2) siny - (63eV/2 + 15eV/4 + 2e 2 g<^) sin 2^ 

27 3 . 
+— e <psin3<p, 
Id 

w 4 = 1134 + 810e 2 + 513e 4 /4 + 135g + 81e 2 g + 57e 4 g/8 + 2q 2 + e 2 q 2 

+ 15r + 18e 2 r + 15e 4 r/8 + s + 3e 2 s + 3e 4 s/8 - (2997e<^ 2 /8 

+1215eV/8 + 225eV 2 /32 + 81eg^ 2 /2 + 45e 3 g<^ 2 /4 + egV/2 

+9enp 2 /2 + 9e 3 V/8 - 27e<//8) cos<^ + (-270e 2 - 531e 4 /8 

-27e 2 g - 59e 4 g/16 - e 2 g 2 /3 - 6e 2 r - 31e 4 r/32 - e 2 s - e 4 s/6 + 

135e V + 45e V/2 + 9e 2 g^ 2 ) cos 2tp + (27e 3 + 711e 5 /256 + 9e 3 g/4 + 

e 3 g 2 /48 + 3e 3 r/16 - e 3 s/8 - 243eV/32) cos3y? - (9e 4 /8 + e 4 g/16 

5 

-e 4 r/32 + e 4 s/120) cos 4^ + e 5 cos5</? + (6237ep/8 + 5e 3 gV/12 

256 

+2835eV/8 + 1455eV/64 + 189eg<^/2 + 135e 3 g^/4 + 3eq 2 <p/2 

+33er<^/2 + 69e 3 r^/8 + 2es<p + 3e 3 sy?/2 - 243e<^ 3 /4 - 135eV/8 

-9eg<^ 3 /2) simp - (675eV/2 + 627eV/8 + 63e 2 g^/2 + 15e*qip/4 

+e 2 q 2 (f/3 + 9e 2 rip/2 + 3e 4 r^/8 - 18eV) sin2y? + (891eV/32 

3 

+135eV/64 + 27e 3 qtp/lQ - 9e 3 r<^/32) sin 3<p - -eV sin 4tp. (90) 

Only the function contains the parameter s. 

The perihelion, as well as, the aphelion of the test particle orbit are 
given by extremals of the function u or w. In (|9*0"D we chose the constants of 
integration to have ip = in the position of the perihelion at any order in e. 
The aphelion is nearby ip = it and the next perihelion closed to ip = 27r. The 
integration constant e is the orbit eccentricity. 

We shall begin the computation of the perihelion shift in the lower 
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order in e. We have 

w = Wq + ewi, (91) 

where Wq and W\ are given in fl9"0|). We verify that, w'(0) = 0. Therefore, the 
perihelion is located the in the angular position (p = 0; the next perihelion is 
in (p = 2tt + L, where L is a small unknown quantity. We have by definition, 
w'(2tt + L) = 0. Expanding the function w'(ip) up to the first order in L in 
the neighborhood of 2n, we find 

w'(2rt + L) - w'{2it) = Lw"{2ti) + 0(L 2 ). (92) 

We obtain the advance of the perihelion up to the first order in e, L = 67re, 
i.e., the usual result for a particle orbiting in the Schwarzschild metric. Now, 
w up to the four order in e is given by the series (0). From the previous 
result we know that the series expansion of L begins with the order e. Thus, 
we will look for a constant L of the form, 

L = Lie + L 2 e 2 + L 3 e 3 + L 4 e 4 , (93) 

solution of the equation 

w f (2ir + L) - «/(2tt) = Lw"(2n) + —w"\2tx) + —w""(2n) + 0{L A ). (94) 

2 6 

Since the fifth derivative is of order e we cut the expansion in the fourth 
derivative . Also, we note that w'(0) = w'(2tc + L) — 0. By identification of 
the coefficients of the different powers of e in the series expansion of the two 
member of Eq. (^) we obtain, 

Li = 6n, 

L 2 = i7r(90 + 15e 2 + 4g), 
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L 3 = vr(1620 + 450e 2 + 120g + 20e 2 g + I2r + 3e 2 r)/4, (95) 
L 4 = tt(379080 + 145800e 2 + 7065e 4 + 38880g + 10800e 2 g + 480g 2 

+80e 2 g 2 + 4896r + 2088e 2 r + 384s + 288e 2 s)/96. (96) 

We see that the discrepancies between the different multipole expansions are 
reflected in L 4 via s. So until the order e 3 we have the same perihelion shift 
for the multipole expansions considered in the present paper. Moreover, 
Eqs. fl9H|) and fl96|) with q = r = s = gives us the perihelion shift of 
an orbit in Schwarzschild geometry up to the fourth order in e. We have 
that the contribution to the perihelion due to deformation of the spherically 
symmetric attraction center appears in the second order in e; a fact that is 
well known in Newtonian theory. 

Now, let us return to the difference equation ([731). For ip and 7 we 
use the solutions up to the five order in u. Therefore we have, £7 = and 
6i/j given by (p4j). The equation ( [73] ) reduces to, 



~ ( - - 2u + 6mu 2 + 0(e)) 5u = 5f3{2l 2 - l)^-p 3 e 4 « 5 . 
2 V p I mr 



In this equation we can replace the solution u, or equivalently w, by his first 
order approximation wq = 1 + ecosy?. We find for the lower order in e, 



esin ip(Su)'- 



3m. „ , 2 

ecosy? (1 + ecosip) 

p 



Su = -SP(2l 2 -1) Al e 4( 1+e cos 
ra A p 



(98) 

By letting e — > in the previous equation we recover the equation for the 



differences of circular orbits, Eq. (|79|) . Since Eq. (Pq) is a linear first order 
differential equation, its solution has a simple integral representation that in 
this case is rather useless due to the fact that the integrals are not elementary. 
The nonhomogeneous part of the solution is proportional to e 4 . Fact that 
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confirms that the first discrepancies between the different expansions are 
proportional to e 4 . 

7 Discussion 

From the study of particular cases, we believe, that are representative of 
different possible series expansions of the metric functions for Weyl solutions 
that have the same Newtonian limit. We found a strong indication that the 
use of the right relativistic multipolar expansion to describe a deformed body 
need to be considered only in the case of a strong gravitational regime like 
bodies orbiting close (a few Schwarzschild radius) around a very compact 
object like a neutron star. For usual planetary motion, in particular for the 
perihelion shift, the effect of having different multipolar expansions with the 
same Newtonian limit can be completely ignored. Also when studying the 
stability of orbits (chaos) of particles moving around deformed bodies, one of 
us [0 found no difference in the trajectories for Newtonian and ERQ defor- 
mations. Recently, the same effect was study on the motion of a gyroscope 

0- 
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Quadrupolar Deviation, Solution I 




-0.45 1 1 1 1 1 1 1 1 1 

23456789 10 

Rs/2m 



Figure 1: We present the deviation for t/> 2 = fpi on ^ ne plane i? = vr/2 for 
a — 1, and different values of the constant a: 0.25 (top curve), 0.5, 0.75, 1.0 
(bottom). 
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Quadrupolar Deviation, Solution II 

I I I I -vvvvr-v^tvvi^vv^-rs;-^^! 




-0.08 1 1 1 1 1 1 1 1 1 

23456789 10 

Rs/2m 

Figure 2: We show the deviation of ip2 = fpi 1 on the plane i? = ir/2 for 
different values of the parameter a: 0.25 (top curve), 0.5, 0.75, 1.0 (bottom). 
The case a = 1 corresponds to the ERQ solution. 
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